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Abstract 

We present a power counting to include Coulomb effects in the three-nucleon system in a low-energy 
pionless effective field theory (EFT). With this power counting, the quartet 5-wave proton-deuteron elastic 
scattering amplitude is calculated. The calculation includes next-to-leading order (NLO) Coulomb effects 
and next-to-next-to-leading order (N LO) strong interaction effects, with an estimated theoretical error of 
~ 7%. The EFT results agree with potential model calculations and phase shift analysis of experimental 
data within the estimated errors. 
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I. INTRODUCTION 



In the last few years there have been a spur of activities involving nuclear effective field the- 
ories (EFT) iQLQ]. Use of EFT in nuclear physics is not a new idea. Chiral perturbation theory 
(%PT) is an EFT that has been phenomenologically successful in the one nucleon sector. The 
next break-through came with Weinberg's work on extending %PT techniques to the two-nucleon 
systems 0,|4|]. Following Weinberg's work, EFTs differing in expansion parameter, dynamical 
content, regularization procedure, etc., were discovered, rediscovered and developed further. The 
initial period was devoted mostly to understanding these EFTs by looking at nucleon-nucleon 
elastic scattering amplitudes S 0, S Hi- From these studies, two alternate formulations of the 
nuclear EFT emerged as relevant for calculating nuclear cross sections at low-energies. The first 
EFT is based on Weinberg's original proposal where one expands the two-nucleon potential in 
perturbation, and then solves the Schrodinger equation with this potential for the nucleon wave 
function. This EFT includes nucleons, pions and photons as dynamical degrees of freedom 10,1a]. 
In the second EFT, one expands the scattering amplitude directly in perturbation, through the use 
of Fe ynm an dia gram s. Here one includes only nucleons and photons as dynamical degrees of free- 
dom [ 1 CM. 1 1 lL 1 1 2l 1 1 3fl . This theory is applicable at momenta p smaller than the pion mass m%, which 
is quite suitable for various nuclear reactions relevant for nuclear astrophysics. In this paper, we 
will work with the pionless EFT as described in Ref. 

Now, we briefly describe the nuclear EFT procedure here. The interested reader should look at 
the comprehensive and up to date review in Ref [2] for details. EFT is a useful tool in the study of 
physical processes with a plethora of clearly separated physical scales. This is generally the case 
in the few-nucleon sector at low-energies, where for example, the deuteron binding momentum 
y ~ 45 MeV is smaller than the pion mass m n ~ 140 MeV, which in turn is smaller than the 
nucleon mass ~ 1000 MeV, etc. EFT provides a natural scheme for separating the short 
distance physics from the long distance effects. At external momenta p ~ m n , pion mass m n sets 
the long distance scale and nucleon mass Mn, heavier meson masses set the short distance scale. 
Currently, we are interested at momenta p ~ y smaller than the pion mass. Thus we construct an 
appropriate low-energy non-relativistic EFT where pion effects are part of the high-energy physics. 
The strong interaction of the nucleons is then described by the most general set of multi-nucleon- 
photon local operators Of, respecting the low-energy symmetries, in an EFT Lagrangian: 

l{x) = £c ; o ; W, (l.i) 

i 

where the effects of the pions and other heavier dynamical particles that were "integrated out" of 
the theory are encoded in the "high-energy" coefficients Q's. The dimensionful couplings C ; 's are 
assumed to depend only on the high-energy scales m % , m p , etc., ~ A and they are determined from 
a fit to experimental data. To make any meaningful prediction with the Lagrangian in Eq. Jl.lt . 
one develops unambiguous power counting rules that determine which operators in Eq. (11.11) are 
most important and which are not. Typically only a few operators are required. Thus one can 
predict a multitude of processes once a few unknown couplings are determined from a few low- 
energy experiments. In addition to this, one also requires power counting rules to estimate loop 
diagrams that describe quantum effects. With these power counting rules one finally expresses all 
physical observables in a perturbative expansion of local operators and loops, where the expansion 
parameter is expected to be Q/A with Q ~ p, y and A ~ m n . The perturbative description of the 
low-energy physics then allows a systematic estimation of errors at any order in the perturbation. 

A word about the high-energy cut-off A is appropriate here. A priori one cannot determine the 
exact value of A in an EFT calculation. This would require complete knowledge of the high-energy 



2 



couplings Q's and all the loop effects, which one does not have. A can be empirically estimated 
from EFT calculations by comparing the contributions from different orders in the perturbation. A 
is process dependent, however, it is found to be ~ m % from various other EFT calculations. 

Recently, effort has been directed towards applications of pionless nuclear EFT, especially in 
the two-nucleon systems involving external currents H13L fl4L ll5L flo. fivl ll8L fl9L 12011 . Progress has 
been made and these calculations have added much to our understanding. In terms of accuracy, 
some of these calculations are at least as good as traditional model calculations. On the other 
hand, the three-body EFT calculations have been so far confined only to the n-d system ifTll EH 
I22L I23L 124 125I1 . In the n-d doublet channel there is a three-nucleon contact operator at leading 
order (LO), whereas in the quartet channel such three-nucleon operators do not contribute up to 
next-to-next-to-leading order (N 2 LO). These calculations reproduce available experimental data 
where applicable within the theoretical errors assumed. 

Application of EFT to study p-d scattering is a natural extension of the work carried out so far. 
An interesting aspect of calculating p-d scattering amplitude would be to study its analyzing power 
A y . This might shed some light on the long standing A y puzzle, p-d scattering will also play an 
important role in EFT calculations of low-energy processes such as pd — > y 3 He, dd —> n 3 He, etc., 
that are important input for primordial light element prediction in big-bang nucleosynthesis [26]. 
Triton beta decay and certain neutrino-deuteron scattering processes receive contributions from the 
same axial current operators with undetermined coupling H15L 12711 . The unknown is one 
of the large sources of uncertainty in the EFT neutrino-deuteron scattering calculations. One could 
in principle determine L\ ^ from triton beta decay data with reasonable accuracy. Understanding 
Coulomb effects in the p-d system will be crucial for the EFT triton beta decay calculation. 

The pionless EFT would be an ideal tool to calculate these low-energy cross sections in a 
model-independent way and to, possibly, reduce the theoretical errors as has been done for the 
np — > dy process H 1 9fl . However, all these require a systematic handling of the Coulomb photons 
in the many-nucleon system as has been done for the purely strong interaction. Elastic p-d scatter- 
ing in the quartet channel provides a unique situation to study Coulomb effects in many-nucleon 
systems. It is complicated enough in the sense that it involves strong interactions and Coulomb ef- 
fects in a few-nucleon system. On the other hand, one does not have to worry about three-nucleon 
forces in the quartet channel as shown in Ref. II llLl2lL 12211 . 

One of the primary goals of this calculation is to establish and understand the EFT power 
counting for the dominant Coulomb corrections at low-momentum. This is crucial for the three- 
nucleon EFT calculations involving more than the n-d system. EFT might play an important role in 
precision calculations of inelastic three-nucleon processes involving external currents. We account 
for the infrared divergent Coulomb contributions and as a first step reproduce in EFT formulation, 
potential model results that have been known for decades. The power counting for the strong 
interaction is well established H3, EH fl3l Eill and it has been successfully applied to calculate 



various two and three-nucleon processes HI ,1111111 111 IIS Mill 1181 19, 20,|2H,|22l|23i 24]. 
First, we recapitulate the strong interaction power counting, ignoring the Coulomb corrections in 
Section HII Our calculations will closely follow the power counting in Ref. |1(1 12, 13, 28j]. Then 
we develop the power counting for the p-d system interacting only through Coulomb photons 
in Subsection IIII Al In Subsection IIIIB1 power counting for the p-d system interacting through 
both the strong and Coulomb interaction is developed. The phase shifts for the quartet 5-wave 
p-d scattering are considered in Section [IVJ We discuss the theoretical and numerical errors in 
the calculation. A comparison with a potential model calculation and phase shift analysis from 
experimental data is also made. Finally, we present our conclusions in Section IVl 
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H. STRONG INTERACTION POWER COUNTING 



The st rong an d Coulomb interactions in the p-d system are described by the low-energy La- 
grangianOMlillii: 
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where "• • • " represents higher dimensional operators with more derivatives. The covariant deriva- 
tive is: 
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and the Pj matrices are used to project on to the 3 5i state, 
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The matrix a, acts on the nucleon spin space and x ; acts on the nucleon isospin space. N is an 
isodoublet field representing the nucleons, andMyv ~ 938.92 MeV is the isospin averaged nucleon 
mass. The auxiliary dibaryon field di has the same quantum numbers as a deuteron and in the 
quartet channel w = — 1. A GauBian integration over the field di in the path integral, followed 
by a field redefinition, reduces Eg. (12.11) to the more familiar nuclear EFT Lagrangian with four- 
nucleon interactions, etc., 112 lL l22l 12911 . The renormalized couplings y and can be determined 
from nucleon-nucleon scattering in the triplet channel 3 5*i II llLl2lLl22l 12911 . 

In the EFT power counting, the expansion parameter is Q/A II 13l 12811 . All physical observables 
are expressed as a perturbation in Q/A. The external momenta p, the deuteron binding momentum 
y and the renormalization scale /u are formally considered 0(Q) and A ~ m n for this low-energy 
EFT. It is assumed in the power counting y 2 ~ 1/A and ~ Q. 

Compared to the LO, we will keep strong interaction corrections up to 0(Q 2 /A 2 ), i.e. N 2 LO. 
Formally, m K /M^ is taken to be o(Q/A), which is numerically consistent. Thus, relativistic 
corrections which typically enter as p 2 / M 2 ,, 'f'/Mjj ~ Q 2 /M^ = Q 2 /m 2 x m 2 /Mfj contribute at 

N 4 LO, and we ignore them here JHElH. 

In this calculation we use dimensional regularization. Some immediate consequences of the 
power counting are, after we integrate over the energy component qo of a loop momentum, con- 
tracting a nucleon propagator: 

1. The loop integration measure / d 3 q scales as 0(Q 3 ). 

2. The nucleon propagator scales as O (M^/Q 2 ). 

3. The dibaryon propagator is —i/<5d at LO. Kinetic energy is ~ Q 2 /Mpj and it contributes at 
next-to-leading order (NLO) and higher. 

Equivalently, every nucleon propagator scales as 0(M^/Q 2 ) and the integration measure 
/ dqod 3 q scales as O [Q 5 / (4%Mn)~\ . We include a factor of 1 / (4tc) with every loop. A closed 
nucleon momentum loop, Fig.[T](a), scales as O [QMn / (4%)]. Thus a two-nucleon loop together 
with a dibaryon propagator, Fig. [T](b), scales as: O [QMn/(4%)] from the loop and a factor of 



4 



o o 

(a) (b) 




FIG. 1 : The single lines represent nucleons, the filled double lines are dibaryons and double line dressed 
dibaryon. The vertex couplings in diagram (a) are 1. All other vertex couplings are y. 

y 2 ~ 1/A from the vertices and a factor of i/o^ ~ l/Q from the dibaryon propagator, which gives 
an overall factor of QM^y 2 / (An<5d) ~ 0{\). Therefore, every diagram can be dressed up by an 
arbitrary number of two-nucleon bubbles with a dibaryon propagator. 

Now, the fully dressed deuteron propagator is given by the sum of diagrams in Fig.[T](c). We 
get, for initial and final spin index {i, j} respectively, 

- -M^ 0- WP V4-AW O -* <2 ' 4) 

This dressed dibaryon propagator represents the deuteron propagator. It is possible to include the 
dibaryon kinetic operator to all orders This allows one to trivially include N 2 LO cor- 

rections due to the effective range p and greatly simplify the calculation. However, the calculated 
p-d scattering amplitude will also include certain N 3 LO and higher order effective range correc- 
tions which should not be included in the strict perturbative sense. Note that one could modify the 
power counting to formally count pA as O (A/ Q) and includ p to all orders in perturbation, as done 
below llioll . We further add that only in the quartet channel (without the three-body force) resum- 



ming the effective range to all orders reproduces the experimental results accurately 1 23. 1241 12511 . 
We find JH: 

1 ij[P0 ' P) l M N f ^p2 /4 _ MnPq _ ie _ _g_ {pQ - p 2 /{4Mn )) 

= iA(p ,p)dij- (2.5) 

The two-nucleon scattering amplitude in the 3 5i channel can now be expressed in terms of the 
deuteron propagator as: 

„(*) = -/.A( P o^ 2 /M«,P = 0) = ^^i-^, (2.6) 

where the cotangent of the 5-wave phase shift 8 can be expressed through the familiar effective 
range expansion: 

/ccot5 = -y+|(/c 2 + y 2 )+w 2 (/c 2 +7 2 ) 2 + --- , (2.7) 
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FIG. 2: Quartet S-wave p-d elastic scattering: the pinball diagrams. The solid lines represent nucleons and 
double lines deuterons. The deuteron-two-nucleon vertex coupling is y. 



where the deuteron binding momentum y= y/MpjB with the binding energy B = 2.224575(9) 
MeV iH. P = 1.765(4) fm JH is the effective range, etc. From Eqs. (1231) and (l2~6l) . we get, 
ignoring shape parameter W2, etc., 

9 871 
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and we also define the deuteron wave function renormalization factor: 



Zd = 
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(2.8) 



(2.9) 



through the LSZ reduction procedure. The amputated amplitudes are multiplied by factors of y/Zj 
for every external deuteron propagator. 

Neglecting Coulomb effects, quartet S-wave p-d scattering only involves neutron exchange di- 
agrams shown in Fig. |2j The tree level amputated diagram involves two factors of y from the 
vertices and a nucleon propagator. Thus it scales as y 2 M^/Q 2 . A n-loop diagram in Fig. would 
include an extra factor of Q 5n / (AnM^) n from the integration measure, a factor of Mjf/Q 4 ' 1 from 
the nucleon propagators, a factor of [471/ ' {y 2 M^Q)\ n from the dressed dibaryon propagators and 
a factor of y 2n from the vertices, for an overall extra factor of O (1). Thus all the diagrams shown 
in Fig. El contribute to the scattering amplitude at LO. In this theory higher order strong interac- 
tions include perturbative corrections to the ratio y 2 /<5d in Eq. ( 12.81) . As mentioned earlier, we 
take effective range corrections into account to all order in perturbation by using the relations in 
Eq. ( 12.81) . From Fig. El we get for the purely strong half off-shell amputated scattering amplitude: 

■criv \ -iM N y 2 
iT s (k,p) 



k-p + k 2 +p 2 -M N E T 



-i4n 



%{k,q) 



d 3 q 

(2tt) 3 _ y+ y/3/4q 2 -M N E T - § (3/Aq 2 -M N E T -f 



(q + p/2) 2 + 3/4p 2 -M N E T 



(2.10) 



where the incoming {outgoing} deuteron carries momentum p{k}, energy p 2 /(4Mn) — 
j 2 /M N {k 2 /(4-M N ) —y 2 /M N -\-z}. Similarly, the nucleon carries momentum — p{— k} and en- 
ergy p 2 1 '(2Mat){& 2 / '(2Mjv) — £}• Hence, the incoming deuteron and nucleon are on-shell and 
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the outgoing deuteron and nucleon are off- shell by an amount £ and — £ respectively. The total 
center-of-mass energy Ej of the p-d system is 3p 2 /(4M^) — y 2 /M^. As in Ref. 12111 . we set 
£ = (k 2 — p 2 ) /Mm and then \p\ = \k\ puts all the external propagators on-shell to give the on-shell 
amplitude. 

IE. COULOMB EFFECTS: p AND Q COUNTING 

Previously, when we neglected Coulomb interactions, it was assumed that the external momen- 
tum p and the deuteron binding momentum y are of similar size, i.e. p ~ y. However, as known 
from non-relativistic quantum mechanics, Coulomb effects enter as QtM^ / p and provide the dom- 
inant contribution at low-energies. Thus in estimating loop effects it is necessary to distinguish 
between the two relevant physical scales p and y. In addition to the expansion parameter Q/A 
we introduce a new expansion parameter ~ p/((xM^). In the present non-relativistic theory there 
is no pair-creation of either dibaryon or nucleon fields. However, a dibaryon field does couple 
strongly to two nucleons. Thus for low-energy p-d scattering, all the diagrams include at most 
one dibaryon field at a given time, which can be put on-shell. There are either one or three nu- 
cleon fields, one of which can be put on-shell, the remaining two propagators being off-shell by 
an amount proportional to the deuteron binding momentum y. Thus every loop integration have 
two dimensionful scales p and y, depending on whether it involves a dressed dibaryon field or 
not. Thus, after integrating over the time component qo and putting one nucleon on-shell, a loop 
integral scales as some power of p or Q depending on weather we pick up the Coulomb correction 
or the strong interaction effects. We explain this in more detail in Appendix lAl 

The previous power counting rules are modified as follows in the presence of Coulomb effects: 

1. The loop integration measure J d 3 q scales as either Q 3 or p 3 . 

2. Every nucleon propagator scales as M^/Q 2 . 

3. The dressed dibaryon propagator scales Q/q 2 . So depending on whether d 3 q ~ Q 3 or d 3 q ~ 
p 3 , the dibaryon propagator scales as l/Q or Q/p 2 . 

4. Photon propagator scales as l/Q 2 or 1/ p 2 , depending on whether d 3 q ~ Q 3 or d 3 q ~ p 3 . 

Rule 3. above is actually not different from the usual power counting where one assumes p ~ 
y ~ Q. When Coulomb photons are involved at low momenta one needs to generalize to the case 
p < y ~ Q. An immediate consequence of these rules is that in a loop with only nucleons, all 
the momenta scale only as Q. These rules become clear when we consider some typical Coulomb 
diagrams. Look at the examples in Appendix lAl as well. 

A. Coulomb Ladder in EFT 

It is easiest to start with diagrams without nucleon exchange, Fig These diagrams repro- 
duce the familiar Coulomb ladder diagrams representing interaction of two charged particles with 
masses equal to M N and M d « 2M^. 

The tree level diagram in Fig. 0(a) is proportional to e 2 / p 2 which is consistent with the power 
counting, (b) is ~ e 2 / p 2 x (y 2 Mjj) / (4-nQ), from the power counting. This is consistent with the 
actual calculation, see Eqs. (IA7I) and (IA8I) in Appendix 1X1 These two diagrams contribute the 
usual tree level Coulomb pieces proportional to l/p 2 . However, Fig. [3] (b) is bigger in the Q/A 
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(a) (b) (c) 





(d) (e) (f) 





(g) (h) 



FIG. 3: Double lines: deuterons; single straight line: nucleon fields; wavy lines: photons. The deuteron- 
two-nucleon vertex coupling is y and the photon-nucleon, -deuteron coupling is e. 

counting. We note that these diagrams with the deuteron wave function renormalization factor 
gives the Coulomb potential in momentum space, if we do a low-energy approximation of the 
nucleon loop integral and keep only the LO contribution from the loop. We have: 

Ana 

z rf x[(a) + (b)] = = y c ( 9 ), (3.1) 

where q is the photon momentum. It is reassuring to recover the familiar result through the EFT 
power counting. 

Fig.|3](c)~ a 2 y 2 Mjj/(pQ^), which looks odd since it scales as I /p. However, as we will show 
later it does not contribute to the Coulomb modified amplitude, in perturbation. With a little effort, 
one concludes that a diagram similar to Fig.|3](c) with n > 3 photon propagators is infrared finite 
and at most scale as y 2 M^/Q 2 x (aM^ / 'Q) n , which is negligible compared to the infrared finite 
contributions from Fig. El and such corrections are ignored. An = 3 photon (attached to a nucleon 
bubble) diagram could contribute as y 2 M^/Q 2 x (aM^ / Q) 3 log(p / Q) < 0.02 for momentum p > 
1 MeV. As we mention later, the numerical procedure used to solve for the scattering amplitude 
does not yield reliable results below momentum p « 20 MeV. Thus log(p/<2) corrections are 
ignored in the present calculation as well. 

From the power counting it follows that dressing any diagram by an extra Coulomb photon 
attached to a nucleon bubble as in Fig.|3](b) contributes a factor of aM^/p. Thus the diagram in 
Fig. EJ(d) ~ e 2 /p 2 (y 2 Mj I )/(4iiQ)aM N /p = (b) x aM N /p, and (e) ~ (a) x aM N /p. Dressing by 
an extra photon attached to just a dibaryon field as in Fig. 0(a) contributes a factor of aM/y /p x 
(4ft<2) / (y 2 Mjf) . Thus these are the effective range corrections to the Coulomb photons attached to 
the nucleon bubble, as can be seen from Fig. 0(f), etc. See Appendix IA1 for more details. Finally, 
dressing by two photons attached to the same nucleon bubble as in Fig.|3](c) contributes a factor 
of ~ a^Mfj/Q 2 w 0.02. Thus Fig. 0(g) is a 2% correction to (b), and (h) is a 2% correction to (c), 
etc. Since we work to only N 2 LO in the purely strong interactions, which have an error of about 
3%, we will ignore these w 2% electromagnetic effects. 
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FIG. 4: The Coulomb scattering amplitude ? c . The solid lines represent nucleons, double lines deuterons, 
wavy lines photons. The deuteron-two-nucleon vertex coupling is y and the photon-nucleon, -deuteron 
coupling is e. 




To summarize, the diagrams in Fig.|3](a) and (b) are iterated to all order and they reproduce the 
Coulomb ladder contribution scaling as 1/ p 2 , 1 / p 3 , etc. We include Fig.|3](c), without iteration, 
which scales as I /p. Iterating any diagram by two photons attached to a nucleon bubble as in 
Fig. |3](c) only modifies the coefficient of the l/p, 1/ p 2 , 1 / p 3 , etc., terms by about 2% and we 
ignore such contributions. The contributions from diagrams with n > 2 photons attached to a 
single nucleon bubble is infrared finite (except n = 3) and negligibly small, and we ignore such 
affects as well. 

We define the amputated Coulomb scattering amplitude by the diagrams in Fig.|4|as: 
.aM 2 y 2 (I \ 1 



i<T c (k,p) = i — ^— I --p 
+i2naM N 



2 Vy V (*-p) 2 

d 3 q T c (k,q) 



(2k) 3 _ y+ ^3/Aq 2 -M N E T - § (3/4q 2 -M N E T - f) 

x (^-p)r J -^' (3 - 2) 

VY / (q-p) 

where the energy-momentum kinematics are the same as in Eq. (12.101) . In Eq. (13.21) . we did not 
include the contribution from Fig.|3](c) since for what we calculate later in Eq. (I4.ll ). contributions 
from Fig.|3](c) in Eq. (13.21) and Eq. (13.31) cancel in perturbation. 



B. Coulomb with Strong Interaction 

Most of what was said about the Coulomb photons also hold here. For example, dressing the 
strong scattering amplitude T s on either side by Coulomb photons, as in Fig.|5](a), enhances T s by 
factors of <xM^/ p. However, there are a couple of differences involving single photon exchange 
diagrams, as in Fig. |5](b) and (c). 

From a naive power counting estimate Fig. |5] (b) ~ (y 2 M^/Q 2 ) x aM/Q which is infrared fi- 
nite. It is as large as the NLO strong interaction corrections to T s so one should include it in the 
calculation. However, a straightforward calculation shows that it is a 7% effect. We do not include 
such contributions in the calculation. Due to the absence of this contribution in the present calcu- 
lation, the theoretical error will be around 7%. The diagram in Fig.|5](c) is a bit more complicated. 
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FIG. 5: The full scattering amplitude TfuII- The solid lines represent nucleons, double lines deuterons, wavy 
lines photons. The deuteron-two-nucleon vertex coupling is y and the photon-nucleon, -deuteron coupling 
is e. 



Naively this particular diagram gets equal sized contributions from the q ~ Q and the q ~ p part 
of the loop integration. Power counting indicates a size ~ Ana/Q 2 . However, a more careful 
analysis shows (c) ~ a/Q 2 log(p 2 /Q 2 ) plus other negligibly small constant pieces. For p > 20 
MeV, neglecting the contribution from (c)< 1% will be consistent with the other approximations. 
We drop contributions from q ~ Q but keep the contributions from q ~ p, for the diagrams in 
(c) and the similar one with a photon attached to a nucleon bubble, for computational ease. This 
approximation turns out to be valid when we compare our results with phase shifts extracted from 
experimental data to within the accuracy assumed. Incidentally, this approximation reduces to 
iterating the strong interaction kernel with the coordinate space Coulomb potential V c = a/r. 
From Fig. |5](d), we get for the amputated scattering amplitude: 



i^FulliKp) 



-iM N y 2 



k-p + k 2 +p 2 -M N E T 2 Vy J (k-p) 2 

TFull(k,q) 



1 



x 



d 3 q 

(2rcF - y+y /3/4q2-M N E T - § (3/4q 2 - M N E T - f ) 
An 



{q + p/2) 2 + 3/4p 2 -M N E T 



2naM N 



VY J(q-P) 



(3.3) 



in the presence of the strong and Coulomb effect, with the approximations mentioned above. The 
energy-momentum kinematics are the same as in Eq. (12.101) . Again we drop contribution from 
Fig.|3](c) for the reason mentioned earlier, following Eq. (I3.2I) . 

Note that diagrams similar to Fig. |5] (c) with more than one photon exchanges are infrared 
divergent. Such contributions are kept, though they might be numerically small for momenta 
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p> 40 MeV. We include all infrared divergent contributions from the Coulomb potential V c = oc/r. 



IV. PHASE SHIFTS 



Predicting the differential cross section for elastic p-d scattering, from Eq. i3.3h . for direct 
comparison with experimental data involves solving a multi-dimensional integral equation. In 
this paper, to reduce the problem to solving a one-dimensional integral equation, we calculate the 
Coulomb subtracted phase shifts (after partial wave projections) instead. However, calculating 
phase shifts has the disadvantage that it depends on precisely the definition used since it involves 
subtracting some Coulomb effects from the full amplitude. To compare our results with available 
phase shift analysis, we use the same subtraction as conventionally defined including resummation 
of the LO Coulomb effects which might not be necessary at momenta above say ~ 40 MeV. 
The conventionally subtracted Coulomb effects in the phase shift analysis are the usual Coulomb 
scattering amplitude with correction for the fact that the charge of the deuteron is not concentrated 
at its center but is bound to the position of the proton in the deuteron H33I1 . 

The Coulomb modified phase shift 8/ is defined as: 



8«(Jk) 

FullW 



5 <0 

} Full 



(*) 



8 (0 



• f 2 l0g 



-/-log 



.2M N k 



' 2M N k (/) 
3n 



(k,k) 



(4.1) 



for every partial wave /. However, it is not possible to project out r ^p\i{k 1 k) and (k, k) for any 
partial wave I because the Coulomb photon propagator is ill defined in the forward direction at any 
momentum k. This problem can be avoided following the proposal by Alt. et. al., see Ref. [34] 
and the references there in for details. The idea is simple: One introduces a photon mass X as a 

regulator and calculates Tpli an d Tc^ using any standard technique for short range interaction. 
Then one numerically reduces the photon mass X until the phase shift 8/ in Eq. (14.11) which depends 



on the difference between Si/ and 8 C 



(/) 



reaches a stable value. 



The full scattering amplitudes with the photon mass X can now be projected onto different 
partial waves and the S-wave amplitude is: 

'■ + p 2 -M N E T 



2M N y 2 l 
kp 2 U 



kp 



2„2 



2pk 
dq 



1 



1 

2 Q 



k 2 +p 2 +X 2 
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Q 



y+^3/4q 2 -ME T - 
M N E T \ 
(IP J 



q 2 + p 2 



§ (3/4<7 2 - 
aM N ( 1 
2% VY 



-M N E T 



q 

t) p 

q 2 +p 2 +X 2 ^ 1 
2qp 



with 



dx — - — 

l x + a 



(4.2) 



(4.3) 
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Similarly, the purely Coulomb 5-wave scattering amplitude with photon mass X is: 



(0)^ aM 2 N y 2 fl k 2 + p 2 + X 2 



?r(k,p;X) = — ^- --pUfi 



/ 
Jo 



2pk \y J 2 V 2k P 
dq 



-Y+ V 3 / 4 ? 2 - me t - 1 (3/4# 2 - M N E T - f) P 

*^(H e (-^)- ^ 



Using Eqs. (14.11 ). ( 14.21 ) and (14.41) . we calculate the 5-wave phase shifts. We find that at any given 
momentum k, a photon mass X in the range fc/10tofc/100 gives a value that is numerically stable to 
about within 1-3% for k ^ 20 MeV, with larger errors for smaller momenta. This does not imply 
that the power counting is invalid below 20 MeV. The sizes of the diagrams are still consistent 
with the power counting estimates. The numerical error is strictly associated with the partial wave 
decomposition. The screening photon mass method shows slow convergence to the X = limit. 
For momentum k below ~ 20 MeV, our numerical routine does not converge in the numerical 
sense 1 . We also have a similar sized theoretical error: higher order strong interaction effects such 
as contribution of shape parameter W2, etc., to the deuteron propagator in Eq. (12.81 ) and Coulomb 
effects ~ a 2 M 2 /Q 2 « 0.02 from diag rams such as those in Fig. 13(g). Therefore, we present 
phase shifts only for momenta k > 20 MeV. The results are shown in Fig. where we present 
the phase shift for n-d scattering as well. For comparison, we also show results obtained using 
potential model AVI 8 |36[ 37] and phase shift analysis obtained from experimental data lHH. We 



note that numerical comparison of the purely Coulomb phase shift calculated in EFT and potential 
model W34I1 agree to within 3% for a given photon screening mass X. 

One can see some expected pattern from Fig. |HJ Coulomb effects are important at low-energies, 
accounting for as much as ~ 40% difference between the n-d and p-d phase shifts at momenta 
k ~ 20 MeV. The 40% effect at k ~ 20 is in the Coulomb subtracted phase shift after removing 
similar sized contributions from the total cross section. It is reasonable to assume that at these low 
momenta the perturbation has at least a significantly better rate of convergence after the resumma- 
tion. It would seem that resummation is unnecessary at higher momenta k> 60 MeV where the 
Coulomb effects in the subtracted phase shift is as low as 6%. However, Coulomb effects would 
still be large in the forward direction in the total cross section. It is also evident that p-d phase 
shift, calculated in EFT, potential model and TUNL phase shift analysis, is consistently larger than 
the n-d phase shifts. 

Note that only certain quantities, e.g. conventionally defined 8' = 8^,, — op , have a sensible 



X — ► limit [34]. In order to ensure a sensible X — > limit and to make a meaningful comparison 
with available partial wave analysis, we use the conventional definition in calculation 8^ including 
resumming Coulomb photons which might not be necessary at large momenta. 

The EFT result for p-d phase shift agrees with the potential model result and the TUNL phase 
shift analysis values, within the estimated ~ 7% theoretical and ~ 3% numerical error. On closer 
inspection, it appears that for momenta p ^ 60 MeV, central value of EFT result is consistently 
smaller than the TUNL values, whereas for very large momentum k ~ 100 MeV it is slightly larger. 
On the other hand the potential model results are consistently larger by a similar amount. 



1 Incidentally, the authors in Ref. Il35ll who also follow the screening photon mass procedure find significant numerical 
errors below momentum k ~ 12 MeV. 
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FIG. 6: Quartet S-wave N-D phase shifts 8 in degrees versus center-of-mass momentum k in MeV. Solid 
curve: EFT p-d result; dashed curve: EFT n-d result; diamonds: potential model results from Ref. ^3(h 
below deuteron breakup and from Ref. %3^1 above breakup at momentum k = 52.77 MeV; dots with error 
bars: TUNL phase shift analysis 

Mi. 



V. SUMMARY 

In conclusion, we developed a power counting for the Coulomb effects in the quartet channel for 
p-d systems in a low-energy pionless EFT. The power counting reproduces p-d scattering results 
that have been known in potential models for decades. However, this power counting will be 
crucial for future EFT three-nucleon calculations beyond the p-d elastic scattering process, where 
precision potential model results might not be available or well known. The quartet channel S- 
wave elastic p-d scattering phase shift was calculated both below and above the deuteron breakup 
threshold. Calculations were performed up to NLO in the Coulomb corrections and N 2 LO in the 
strong interactions, with an estimated theoretical error of ~ 7% and ~ 3% respectively. An error of 
about 2% was found in the numerical evaluation of the scattering amplitude. Within the estimated 
error, the EFT results agree with both potential model calculations and phase shift analysis of 
experimental data. 

The LO and NLO Coulomb ladder contribution is equivalent to contributions from the coordi- 
nate space Coulomb potential V c = a/r. The large 7% error is primarily from the single photon 
diagrams, similar to that in Fig.|5tb) dressed with the LO scattering amplitude on the external legs. 
These diagrams constitute the largest higher order corrections to the Coulomb ladder. There have 
been some discussion in the literature about Coulomb polarization effects. These effects have been 
found to be negligible for the present calculation [35] and they are not included. 

The power counting for the Coulomb effects developed here could be applied to higher partial 
waves in both quartet and doublet channel for p-d scattering. As mentioned in the introduction, 
one can then go on to study the so-called p-d A y puzzle, and low-energy many-body processes 
such as pd — > y 3 He, dd — > n 3 He, etc., that might be important for placing stricter bound on 
primordial light element abundances in big-bang nucleosynthesis calculations tl26fl . Triton beta 
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decay is another important low-energy process where three-body Coulomb effects are important. 
This process is related to neutrino-deuteron scattering and might have significance to neutrino 
physics at SNO. 
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APPENDIX A: THE DEUTERON PROPAGATOR WITH COULOMB PHOTON 



To understand the power counting rules of Section Unl it is most convenient to start with the 
leading order deuteron propagator, without the range correction: 

>a( ^ ■ 4n 1 



M N y 2 -y+ v /p 2 /4-M N p -ie 



471 y+ y / p 2 /4-M NP0 



Mn)> 2 p 2 /4 - M N p Q - y 2 - ie ' 

In elastic p-d scattering, every Feynman diagram with a deuteron propagator involves only a 
single nucleon propagator at equal times, see Figs. I2I3I4I and The deuteron propagator mo- 
mentum in any diagram can be written such that (po = Ed + qo,P — q) with the corresponding 
nucleon line carrying energy-momentum (E^ — qo, —q) where Ed = p 2 / (4M#) — y 2 /M^ is the in- 
coming/outgoing deuteron energy and Epj = p 2 / (2Afjy) is the nucleon incoming/outgoing energy 
with incoming/outgoing momenta p, for a generic loop momentum (qo,q). We get: 



4n Y+a/^-A^o+Y 2 



iA(E d + q 0j q) = -i— 5 -|— 2 , (A2) 

M N y l 3_Pt_ MNqQ _ iE 

where the factors of y 2 cancel in the deuteron pole. Carrying out the qo integral, as in Eq. (IA5L 
and putting the nucleon on-shell with qo — E^ — q 2 / (2M N ) gives 



r_ = 4k Y+\/3^+Y 2 
IM n l M N y 2 

Including the effective range p is straightforward, and we get 



iA{Ed+EN ~^ q) = "W ^EZ ■ (A3) 

4 



q 2 An Y+V3^+Y 2 1 

iA(E d + E N --?—,q) = -i—^j- V . \ . . (A4) 

™* M N y 2 3 <t=£ 1 p/ J W _ P 2 — \ 



~ 1-5 Y+V3^+Y 
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We will ignore the effective range corrections ~ py ~ Q/A in the following discussion for simplic- 
ity. For small momenta p <C y ~ Q, when q ~ p (d 3 q ~ p 3 ) we get A ~ (4k) / (M^y 2 ^/ p 2 ~ Q/p 2 
whereas when q ~ y (d 3 q ~ y 3 ) we get A ~ (4tt) / (M^v 2 y) ~ 1/2. This scaling of the deuteron 
propagator is far from obvious. Even though one might naively expect A ~ l/Q from dimensional 
analysis in Eqs. (12.41 ) and (12.81 ). this is invalidated since the deuteron can be put on shell which 
exactly cancels the factors of y 2 from the deuteron pole. 

Now, even though the contribution from the deuteron pole is enhanced for q ~ p, its contribu- 
tion is typically suppressed by d 3 q ~ p 3 except when the loop integral also involves a Coulomb 
photon propagator which scales as l/q 2 ~ l/p 2 - This explains why there are no l/p infrared 
enhancements in the diagrams without photons in Fig.|2j where the dominant contribution comes 
from loop momentum q ~ y ~ Q with A ~ l/Q. On the other hand diagrams with Coulomb pho- 
tons get infrared enhancements from q ~ p, and A ~ Q/p 2 - For example, from Fig. [21(f) we get 
the contribution (ignoring factors of z',2, K, etc.,): 

dqo d 3 q 1 a a 



271 W -qo + E N -^(q-p) 2 (q-ky 

d 3 q k . q 2 , a a 

*-A(E d + E n - v 1 



(271)3 * " 2MN^(q-p) 2 (q-k) 2 
a 2 f d 3 q (y+v/3^ 2 /4-3p 2 /4 + y 2 ) 1 



My 2 J (2k) 3 3q 2 /A-3p 2 /A (q-p) 2 (q 2 -k) 



(A5) 



where the deuteron energy E c i = p 2 / (4Mpj) — y 2 /Mjy and the nucleon energy E^ = p 2 / (2M^) with 
momenta p, k defined as in Eq. (|2.10b . This integral is dominated by momenta q ~ \p\ — \k\ and 
the contribution is, without wave function renormalization, 

a aQ a aM N Q 

~ — n 771 — = 9 X - . (A6) 

p L My L p p L p Myy 1 

Diagrams without a deuteron propagator involve three nucleons. After carrying out the integral 
over the energy component qo and putting a nucleon on-shell, one is left with two nucleons that 
are always regulated in the infra red by the deuteron binding momentum y, as these involve two 
nucleons that carry the energy Ej of the "parent" deuteron. For example, in Fig.|3](b) we have 



)' 2 



dqo d 3 q 11 la 



2k (2k) 3 qn i 2 a - I Ea (g+p)2 a,i I E, {q+k)2 (P~ k)2 
ay 2 M 2 f d 3 q 1 1 



p-k) 2 J (271) 3 (q+P) 2 I Q 2 P 2 I V 2 (q+k) 2 , q 2 P 2 , .,2 
2^2 2" 1 "' 2" r 2 2" 1 "' 



ay 2 A/ 2 



P 2 e 



(A7) 



where the nucleon propagators get regulated in the infrared by the deuteron binding momentum 
y ~ Q with \p\ = \k\. After wave function renormalization Zj ~ Q/ (M^y 2 ), it gives a contribution 
~ a/p 2 . Exact evaluation of the integral gives 



1 

arctan 



2K( P -k) 2 



(p-k) 2 
2 V 4y 2 - /? 2 



g^+°(A (A8) 
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which is consistent with the power counting estimate (including all the numerical factors). 

To summarize, all the nucleon propagators scale as M^/Q 2 , the deuteron propagator scales as 
Q/q 2 and the Coulomb photon propagators scale as l/q 2 . The loop momentum scales as either 
q ~ p <C Q or as q ~ Q, and correspondingly the deuteron propagator scales as l/Q or Q/p 2 . 
In estimating the dominant contributions from a loop integral for q ~ p, one has to consider the 
suppression factors from d 3 q ~ p 3 along with infrared enhancements from the deuteron Q/p 2 and 
the photons 1 / p 2 . 
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